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We present an analytic theory of the spin-resolved pair distribution functions g arT '(r) and the 
ground-state energy of an electron gas with an arbitrary degree of spin polarization. We first 
use the Hohenberg-Kohn variational principle and the von Weizsacker-Herring ideal kinetic energy 
functional to derive a zero-energy scattering Schrodinger equation for -J ' g aa i{r). The solution of 
this equation is implemented within a Fermi-hypernetted-chain approximation which embodies the 
Hartree-Fock limit and is shown to satisfy an important set of sum rules. We present numerical 
results for the ground-state energy at selected values of the spin polarization and for g aa i (r) in both 
a paramagnetic and a fully spin-polarized electron gas, in comparison with the available data from 
Quantum Monte Carlo studies over a wide range of electron density. 
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03 ■ I. INTRODUCTION 

The homogeneous electron gas (EG) of fermions interacting by the e 2 /r law and moving in a uniform neutralizing 
background of positive charge has been for many decades the basic reference system for calculations of electronic 
structure in condensed matter physicsi^ The extent to which exchange and correlations compete with each other on 
varying the electron density has been a central issue going back to the work of Bloch on spontaneous spin polarization 
due to exchange and to that of Wigner and Seitz on cohesion in simple metals. The many-body effective interactions 
I ' acquire a magnetic component arising from local deviations of electron density around any given electron, with the 
, exchange term being determined by the electrons carrying the same spin and the correlation term being to a large 
t-H ' extent determined by the electrons of opposite spin. The role of spin fluctuations and their correlations thus are a 
subject of continuing interest, on which theoretical progress starting from the basic Coulomb Hamiltonian has been 
C^J painfully slow. 

There is a broad consensus on the increasing importance of spin polarization at strong coupling and of correlations 
in lowered dimensionality. This fact has been repeatedly emphasized by Quantum Monte Carlo (QMC) numerical 
studies, starting from the early work of Ceperley and Alder£ In particular, in the three-dimensional (3D) electron 
system these studies have revealed a continuous transition from a paramagnetic to a fully spin-polarized ground state 
occurring with increasing coupling strength, before a first-order quantum phase transition to a ferromagnetic Wigner 
crystal occurs. 4 5 While the 31? system appears to be most suitable at present for improvements and critical tests of 
I ■ analytic theories, it may be mentioned that studies of spin motions in structures of reduced dimensionality are being 
greatly stimulated by expectations of developments in spin-based electronics and in quantum computing^ 

The main purpose of the present work is the formulation of a practicable analytic theory for spin correlations in 
the 31? EG and the comparison of its numerical predictions with the available QMC data. The fluid that we consider 
consists of two spin components at fixed average densities n a with a =f or J,, and the spin correlations are described 
by the pair distribution functions g a(T i{r). These are defined so that the quantity Amr^n^i g aa i (r)dr gives the average 
number of electrons with spin a' lying within a spherical shell of radius r and thickness dr centered on an electron 
with spin a. The general strategy is to set up a Schrodinger-like differential equation for \/ g a a' ( r ) with the help of the 
Hohenberg-Kohn variational principle^ and to solve it by a Fermi-hypernetted-chain type of approximation tailored 
to embody the Hartree-Fock limit as well as a set of sum rules for the two-component Fermi fluid. 

It may be mentioned at this point that there has quite recently been a renewed interest in the study of the spin- 
averaged pair distribution function g(r) in paramagnetic EG models within a two-body scattering approach first 
proposed by Overhauseri^ In brief, a Schrodinger equation is set up to evaluate particle-pair wave functions from 
which g(r) can be built through sums over angular momentum and spin. Various approximations have been examined 
for embodying the many-body effects into the effective scattering potentials^ In an alternative approach, which is 
fully developed in the present work, we have instead sought to use a density-functional framework for a self-consistent 
evaluation of g(r) and of the effective scattering potential. 9 

The contents of the paper are briefly as follows. Section II presents a formally exact zero-energy scattering equation 
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for the partial distribution functions g aa i (r) and introduces the approximations that we propose for the evaluation 
of the three main contributions to the scattering potential. Section III shows that these approximations satisfy the 
plasmon sum rule, the charge neutrality condition, and Kimball's cusp condition. Our numerical results for the 3D 
EG are presented in Section IV. Finally, Section V summarizes our main conclusions and gives suggestions for further 
work. 



II. THEORY 



We consider an inhomogeneous 3D fluid of electrons consisting of the two spin species with densities n a (r) in 
the presence of external potentials V° xt (r) and of a uniform neutralizing background. From the Hohenberg-Kohn 
theorem^ the ground-state energy functional of the fluid can be written as 

£ gs [{n CT (r)}] = T s [{n a (r)}] + £ f d 3 r ^ cxt (r)An CT (r) + £„[{n CT (r)}] + E KC [{n a (r)}] (1) 

where An CT (r) = n CT (r) — n a are the deviations of the spin densities from their average values, T s and E xc are the ideal 
kinetic energy and exchange-correlation energy functionals, and E H is the Hartree term given by 

E H [{n a (r)}} = \j^J <?r J d 3 r' v(\r - r'|) An CT (r)Arv (r') (2) 

with v(\r — r'|) = e 2 /|r — r'|. The presence of a neutralizing background has been taken into account in these equations. 
Following Herringifi we decompose the kinetic energy functional into the sum of two terms, 



where the first term is the von Weizsacker "surface" kinetic energyJA and Te is itself defined by Eq. (3). 

This general formalism is adapted to the derivation of differential equations for the spin-resolved pair distribution 
functions g atT > (r) in the homogeneous EG by viewing the quantity n a \g aa ' (r) — 1] as the distortion that an electron 
in the EG, located at position r — with spin a' , induces in the density profiles of the electrons with spin ail The 
appropriate ground-state energy functional for this problem is obtained from Eqs. (l)-(3) by the formal replacements 



m — > /i 

Aext 



V CT ext (r) ->u(r) 
n a (v) -> n a gaa'(r) 



(4) 



where /i = m/2 is the reduced mass of an electron pair. The Euler-Lagrange equations for the spin- resolved pair 
functions can now be obtained from the variational principle of Hohenberg and Kohn£ using the von Weizsacker- 
Herring ideal kinetic energy functional as shown in Eq. (3). With the zero of energy taken at the chemical potential, 
the formally exact differential equation for g aa i (r) reads 



Vl + v{r) + v% a (r) + V™ (r) 



(r) = . 



Here, the "Pauli potential" Wp CT (r) is defined by 



v7 (r) = 



5Te[{n a (r)}} 



5n a (r) 



n CT (r)=n„ g CT<T / (r) 



and (r) is the "excess" Kohn-Sham potential, which is given by 



VZ° (r) - Mr) + <° (r) = n d 3 r'v(\r - r'\)[g(r') - 1] + 



6E* 



>(*)] 



8n acT > (r) 



(5) 



(6) 



(7) 



n a (r)=n a g aa i{r) 



with n = ti t +n l and g(r) = J2a,cr'( n <? n <r' l n2 )9<j<j>{r). 
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Equations (5)-(7) show that a zero-energy scattering theory approach to pair correlations in a quantum electron 
fluid has a sound theoretical justification within the framework of the Hohenberg-Kohn variational principle. Of 
course, the functional dependence of Tg and of £ xc on the inhomogeneous electron densities is not known and we shall 
need to resort to appropriate calculational schemes and approximations for the potentials introduced in Eqs. (6) and 
(7). 

A possible way to handle the kinetic energy term is to pass to a Kohn-Sham scheme, by expanding the pair 
distribution functions into Kohn-Sham two-particle scattering orbitals. This was done in a series of previous studies^ 
as already mentioned in Section I, but the model scattering potentials used there fail to satisfy the plasmon sum rule 
and therefore do not account correctly for the behavior of g CT(T /(r) at large r. In the following we shall instead take 
advantage of the formal similarity between Eq. (5) and the Euler-Lagrange equation which is adopted in the so-called 
Fermi hypernetted-chain approximation (FHNC). This is derived by a Jastrow-Feenberg variational Ansatz on the 
many-body wave functioniiiii^ and was used to treat two-component Fermi fluids in the case of electron-hole liquids 
and of liquid metallic hydrogen^ Such an approach will allow us to build into the theory some important sum rules 
and limiting behaviors. 



A. The Pauli potential 



In the FHNC context the Pauli potential can be chosen so as to ensure that the Hartree-Fock limit is correctly 
embodied into the theory. An argument can be given to show that this term in Eq. (5) becomes dominant at weak 
coupling by considering how the various functionals scale under scaling of all lengths by a factor A. It was shown 
by HerringiS that Tg scales by a factor A~ 2 , whereas Levy and Perdewii have shown that both the Hartree and 
the exchange energy functionals scale by a factor A -1 . Only the inequality E c [{n a (r)}} > A -1 E c [{n^\r)}} could be 
proved^ for the correlation energy functional under scaling, but it seems reasonable to expect that on taking A = r s 
the correlation functional would also become relatively negligible in the limit r s — > 0. Here r s is the usual coupling 
strength parameter for the EG, which in 3D is related to the average electron density n by r s as — (47rn/3) -1 / 3 with 
as the Bohr radius. 

From the above argument it follows that Eq. (5) implies 



h 2 vgyggMO 



(8) 



in the weak coupling limit r s 
approximation (HF), namely 



0. In Eq. (8) g™,(r) are the spin-resolved pair functions in the Hartree-Fock 



<C(r) = l-9 [h{k Fa r)/{k Fa r)Y 
9%{r) = 1 



(9) 

n(cr)C) 1/3 with 



where a — —a, ji(x) = (sinx — xcosx)/x is a spherical Bessel function, and kp a = &f(1 + sg 
kp = (Sn 2 !!,) 1 ^ 3 the usual Fermi wave number and £ = |n| — n\\jn the degree of spin polarization. 

Although the expression for the Pauli potential in Eq. (8) is correct only for a weakly coupled Fermi fluid, we 
shall assume in the following that it can yield useful results in our self-consistent calculations of the pair distribution 
functions with increasing coupling strength. This assumption will have to be tested through quantitative comparisons 
of our numerical results with the available QMC data. As a broad qualitative argument in support of this assumption 
we may remark that the role of the statistics is expected to weaken with increasing coupling strength at least in the 
spin-averaged pair functioni&iS and that only the von Weizsacker term contributes to the kinetic energy functional 
in the case of a charged Bose fluids 



B. The excess Kohn-Sham potential 

The FHNC expresses the potential V^ 7 (r) in Eq. (5), which is the sum of the Hartree and of the exchange- 
correlation potential, as the sum of two effective pair interactions ii&i^iSiiii 

VZ a '{r) = W^\r) + <WWT(r) • (10) 

The first term on the left-hand side of Eq. (10) descends from the two-body correlation functions u aa i(r) in the 
Jastrow-Feenberg wave function and is therefore formally the same as for a binary boson mixture. The second term 
is instead due to the antisymmetry of the fermion many-body wave function. 
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As shown by Chakrabortyi^ in treating a binary fermion mixture, the HNC closure yields for u aa > (?*) the expression 

Uaa'(r) = lnflw/(r) - [g aa -{r) - 1] + <w(r) (11) 

where c CTCT / (r) are the direct correlation functions, which are related to g a(7 ' (r) by the Ornstein-Zernike relations. We 
introduce at this point the partial structure factors S aa ' (q) of the binary mixture, which in essence arc the Fourier 
transforms of g aa i (r) : 

Saa'(q) = <W + yjn a n a , J d 3 r \g aa >{r) - 1] exp(-iq ■ r) . (12) 

We also introduce the Fourier transform of (r) (W£ a (q) , say). Minimization of the ground state energy against 
arbitrary variations of g aa <(r) yields with the help of Eq. (11) the expression 

W^'(q) = --^[S aa> {q) - <W] - V aa ,(q) (13) 



where e q = h q /(2m) are the single-particle kinetic energies and the functions V acr i(q) are given by 

V.Aq) = {-1 + [SUl) + S^(q)]/A 2 (q)} 
V as (q) = - fg_ g gff (g) [S aa (q) + S^(q)} /A 2 (q) 



(14) 



with 

A(q) = S air (q)S 99 {q)-& aS {q). (15) 

Equations (13)-(15) show how the effective boson-like interactions (r) in Eq. (10) are related in Fourier transform 
to the spin-resolved pair distribution functions. 

Turning to the second term on the left-hand side of Eq. (10), the effective pair potential W° a (r) has a very 
complicated expression within the FHNCiiiiiiiii However, in dealing with a one-component electron fluid Kallio 
and Piilo^S have proposed a simple and effective way to account for this consequence of the antisymmetry of the 
fermion wave function. Their argument is immediately generalized to our two-component Fermi fluid, and leads to 
the requirement that in Fourier transform this term should cancel the effective boson- like interaction W£ a (q) for 
parallel-spin electrons at low coupling. That is, 



Wr(q) = Km Wr(q) = £-[1 + 2 S£(g)] 



(16) 



Here, S^(q) is the Hartree-Fock structure factor, which is given by 



chf/ x _ / f (<z/fcpv) - jeiq/kFa) 3 for q < 2k Fa . . 

O tr aW-^ 1 forq>2fc Fff . [ <> 

It is evident that the insertion of Eqs. (8)-(17) into Eq. (5) allows a self-consistent calculation of the spin-resolved 
pair distribution functions and of the effective electron-electron interactions. Before proceeding to the numerical 
solution of this problem, we examine how the approximate theory presented above fares in regard to some exact 
properties of pair correlations. 

III. SUM RULES AND LIMITING BEHAVIORS 

In this Section we show that the pair functions obtained from the theory presented in Section II satisfy three exact 
properties. These are the charge neutrality condition, the plasmon sum rule, and Kimball's cusp condition. 

The asymptotic behavior of the effective potential V aa '(r) is first obtained from Eq. (5), which can be rewritten in 
the form 

V aa .{r) = 9aa'(r)[v(r) + HT>) +<">)] + fcw(r) - W'W + ^ ^V^F)]" ■ (18) 
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A careful analysis of this equation shows that V aa '(r) — » v{r) + 8 aa iW"' T {r) for r — > oo, and hence in Fourier transform 
we have from Eq. (16) 



V aa {q) -^v q + 8ir 2 h 2 /(3mk F a) 



(19) 



for q — > 0, with v 9 = 4ire 2 /q 2 in the 3D EG. The corresponding asymptotic behavior of the structure factors is 
obtained from Eqs. (14), which can be inverted to yield 



S<j<j(q) = 



A** A 2 - 



ALA 3 



2+[A a 



A SS \A 



(20) 



where A (7(7 '(q) = (W + 2sgn(aa')^/n a n a > 
find 



^2 + [A aa + A SS ]A 
Vaa'{q)/e q and A(q) = (A n A u - A 2 n )~ 1 / 2 . From Eqs. (19) and (20) we 



(21) 



for q-^Q, with w p i = {Anne 2 /m) 1 / 2 , g aa {Q = [1 + sgn(<r)C], g aS {Q = V 1 ~ C 2 , and 

:! /2[l-sgn(a)C]/[l + sgn(a)C] 



/-(C) - 



U(0 = 



(i + C)- 1/3 + (i-C)- 1/3 



(l+o-vs + a-o- 1 / 3 



(22) 



These limiting expressions were earlier known to hold at £ = within the Random Phase approximation (see e.g. 
Gori-Giorgi et alm^ and references given therein). The form of Eq. (21) immediately ensures that the charge neutrality 
condition is satisfied. This reads 

-<W (23) 



(24) 



y/n a n a > J d r [sw(r) - 1] = -<W 

or 5 ffo -/(g -> 0) = from Eq. (12). 

The form of Eq. (21) also ensures that the plasmon sum rule is satisfied. This reads 

lim S(q) = e q /(hw p i) 

q^0 



where S(q) = YL a , a ' [\J n <? n <y' /n)S a „i (q) . 

The cusp condition as first demonstrated by Kimball 22 relates the logarithmic derivative of g^ (r) in the origin to the 
Bohr radius and is satisfied if the scattering potential between pairs of electrons with antiparallel spins reduces to the 
bare Coulomb potential for r — > 0. This is ensured in our approach by the properties v%° (r — > 0) = W" a (r — > 0) = 
and W^ a (q — > oo) oc q~ 6 , with the latter form being equivalent to the asymptotic behavior Sn(q — > oo) oc q~ 4 . 

Finally, the Pauli principle requirement g CTCT (0) = is ensured in our approach by the behavior of the Pauli potential 
in the origin, v° a (r —> 0) — ► 2 S 2 /(mr 2 ). This behavior determines the power-law exponent in the distribution function 
of parallel-spin electrons, g a a{f — > 0) oc r 2 , as can be proved directly from Eq. (5). 



IV. NUMERICAL RESULTS 



We turn to a presentation of our numerical results, which are obtained by solving Eq. (5) with the following self- 
consistency cycle. We start with the trial choice g a a'{r) — g^,{r) and W^ a (r) = 0, and find the effective potentials 
V aa ,(r) by means of Eq. (18) and hence the structure factors S aa *{q) via Eq. (20). At this point we can calculate 
new values for g a cr'{r) and for W^ a (r) by taking Fourier transforms and using Eq. (13). This procedure is repeated 
until self-consistency is achieved. The computational time typically needed to obtain g aa i (r) at each value of r s and 
£ is a few minuted on a PC with a Pentium IV/1.4 GHz processor. 
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We have calculated in this way the spin-resolved pair distribution functions of a 3D electron gas for r s up to 100. 
The main results of our work are shown in Figures 1-8. 

In Figure 1 we show that our results for g(r) = [g^(r) + gn(r)}/2 in the paramagnetic EG at r s = 1, 5, 10 and 20 
are in excellent agreement with the QMC data of Ortiz et al.m In the same range our results for g(r) are in excellent 
agreement with those of Kallio and Piilo. 20 In Figures 2 and 3 we show that for the same cases our results for the 
spin-resolved pair functions <?tt( r ) an d 9tl( r ) are a l so m excellent agreement with the QMC data of Ortiz et al.. A 
To the best of our knowledge theoretical results of similar quality have not been reported in the literature from an 
approach which is free of input and/or fitting parameters. Figure 4 shows a comparison of the present results at 
r s = 10 with those previously obtained by us& with a self-consistent Hartree approximation. 

In Figure 5 we compare our results for the paramagnetic EG at r s — 40 and 100 with the QMC data of Ortiz et al~ 
and with those obtained more recently by Zong et al.A As discussed by the latter authors, these two QMC approaches 
are essentially the same but involve some technical differences. In particular the study reported by Zong et al& 
includes backflow and three-body terms in the wave function, and uses a novel numerical technique termed "twisted 
averaged boundary conditions" which allows a sizeable reduction of finite-size errors^ 3 - The reader is referred to Ref. 
I^for details. These differences are not expected to change significantly the results of Ortiz et al. 4 for r s below 20, but, 
as is clear from Figure 5, the differences between the two QMC studies become significant at larger r s . Our results 
lie somewhere in between the two sets of QMC data at large r s . From the theoretical point of view it is important to 
remark that in FHNC-type calculations at strong coupling the induced W^ a (q) interactions should be corrected by 
the addition of three-body correlations and elementary-diagrams (or "bridge functions") contributions i^i^iik 

Further comparisons with QMC data on pair distribution functions can be made for the fully spin-polarized EG 
(C = 1). This is done in Figure 6 at r s — 1,4,40 and 100. No QMC data are available in the intermediate range 
4 < r s < 40, but the good agreement with the QMC data of Ortiz et alA in the low-r s range is evident. The same 
inaccuracies in the theory that we have exposed in the paramagnetic EG at strong coupling are also found at Q = 1. 

Finally, in Figure 7 we report the spin-dependent effective electron-electron interaction V°g (r) = v(r) + W-^ (T (r) + 
S (T(T >Wg' J (r) and its Fourier transform V°£ (q), as it emerges from our self-consistent calculations on a 3D EG both in 
the paramagnetic and in the ferromagnetic state. The attractive part of the parallel-spin effective interaction deepens 
at increased £ as is physically expected. The Pauli potential and its dependence on £ are shown in Figure 8. 



The ground-state energy e g (per electron) of the EG at each value of r s and C can be calculated by means of an 
integration over the coupling constant A, 



directly from the structure factor S\(q) calculated with interaction v q = 47re 2 A/g 2 . Here Sq = 3[(1 + C) 5 ^ 3 + (1 — 
£) 5 / 3 ]/(10a 2 r 2 ) Ryd is the ideal-gas kinetic energy with a = (97r/4) -1 / 3 , and S(q) has been defined immediately 
below Eq. (24). In fact, the integration over A is carried out by integration over r s . The second term in Eq. (25) is 
the exchange-correlation (xc) energy, from which the correlation energy is obtained by subtracting the expression for 
the exchange energy e x = -3[(1 + C) 4/3 + (1 - C) 4/3 ]/(47rar s ) Ryd. 

We have calculated the ground-state energy for ( = 0, 0.333, 0.667 and 1 over the range 1 < r s < 50. The results 
are reported in Table I in comparison with QMC data from Ceperley and Alder ji Ortiz et alA and Zong et al.A The 
Table also includes other theoretical results obtained in the self-consistent dielectric theory of Singwi et alw^ and in 
the modified-convolution approximation of Tanaka and Ichimaru^ 

It is seen from Table I that the present theoretical approach yields fairly accurate values of the ground-state energy 
even at r s ~ 50, even though the details of the pair distribution function are becoming quantitatively inaccurate as 
we have already seen in Figures 5 and 6. Much higher accuracy is needed for theoretical predictions on the quantum 
phase transition from the paramagnetic phase to the fully spin-polarized phase, since from the QMC study of Zong et 
al~ the difference in energy between these two phases at r s = 100 is only about 4 x 10 _6 Ryd in favor of the latter. A 
continuous transition between these two states seems to start for r s ~ 50-70 according to the QMC data. In contrast, 
the theoretical approaches of Singwi et alS^ and of Tanaka and Ichimaru 2 ^ suggest that the magnetically ordered 
state may be the favored one already at r s ~ 40-50. 



A. Ground-state energy 




(25) 
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V. SUMMARY AND DISCUSSION 



In summary, we have presented in this work a new theoretical study of the spin-resolved pair distribution functions 
and of the ground-state energy of the 3-D electron gas. Our approach yields numerical results of good quality in the 
regime of weak and intermediate coupling strength. 

Improvements of the theory will be necessary for a quantitative study of the magnetic phase diagram and for an 
extension to electron fluids of lower dimensionality. As we have already commented in the main text, the kinetic 
energy functional may be treated more accurately by recourse to an analysis of the pair functions into Kohn-Sham 
pair orbitals. 8-9 Preliminary calculations indicate that an appreciable improvement can be obtained in this way in the 
values of the ground-state energy at small to intermediate coupling over the values reported in Table I. At larger values 
of the coupling we expect that a sizeable improvement on the FHNC approximation to the effective electron-electron 
interactions can only arise from the inclusion of contributions from the so-called elementary diagrams and from 
three-body correlations Better quantitative accuracy can be expected to be achieved by such means when the self- 
consistency requirements on the theory are extended to include the thermodynamic sum rules on the compressibility 
and spin susceptibility of the electron fluids 



Acknowledgments 

This work was partially supported by MIUR through the PRIN2001 Initiative. We are indebted to Prof. P. Ballone 
and Prof. D. M. Ceperley for providing us with their QMC data reported in the Figures. 



1 D. M. Ceperley, Nature 397, 386 (1999). 

2 For a recent review see Electron Correlations in the Solid State, edited by N. H. March (Imperial College Press, London, 
1999). 

3 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 (1980); B. J. Alder, D. M. Ceperley, and E. L. Pollock, Int. J. 
Quant. Chem. 16, 49 (1982). 

4 G. Ortiz, M. Harris and P. Ballone, Phys. Rev. Lett. 82, 5317 (1999). 

5 F. H. Zong, C. Lin and D. M. Ceperley, Phys. Rev. E 66, 036703 (2002). 

6 For a recent review of semiconductor spintronics see S. Das Sarma, J. Fabian, and I. Zutic, cond-mat/0206540 see also J. 
Fabien, I. Zutic, and S. Das Sarma, Phys. Rev. B 66, 165301 (2002). For reviews on dilute ferromagnetic semiconductors see 
J. Konig, J. Schliemann, T. Jungwirth, and A. H. MacDonald, in Electronic Structure and Magnetism of Complex Materials, 
edited by D. J. Singh and D. A. Papaconstantopoulos (Springer, Berlin, 2003) Chapter 4; see also T. Jungwirth, J. Sinova, 
J. Kucera, and A. H. MacDonald, cond-mat/0210318 

7 R. M. Dreizler and E. K. U. Gross, Density Functional Theory, An Approach to the Quantum Many-Body Problem (Springer, 
Berlin, 1990). 

8 A. W. Overhauser, Can. J. Phys. 73, 683 (1995); P. Gori-Giorgi and J. P. Perdew, Phys. Rev. B 64, 155102 (2001); B. 
Davoudi, M. Polini, R. Asgari, and M. P. Tosi, Phys. Rev. B 66, 075110 (2002); F. Capurro, R. Asgari, B. Davoudi, M. 
Polini, and M. P. Tosi, Z. Naturforsch. 57a, 237 (2002). 

9 B. Davoudi, R. Asgari, M. Polini, and M. P. Tosi, Phys. Rev. B 67, 172503 (2003). A preliminary report of some results 
for both bosons and fermions has appeared in cond-mat/0206456 (unpublished) and in poster i?1.299 presented at the APS 
March Meeting 2003. 

10 C. Herring, Phys. Rev. A 34, 2614 (1986). 

11 C. F. von Weizsacker, Z. Phys. 96, 431 (1935). 

12 The classical analog of this prescription was first considered by J. K. Percus, in The Equilibrium Theory of Classical Fluids, 
edited by H. L. Frisch and J. L. Lebowitz (Benjamin, New York, 1964), p. 11-83. 

13 L. J. Lantto and P. J. Siemens, Nuclear Phys. A 317, 55 (1979); L. J. Lantto, Phys. Rev. B 36, 5160 (1987). 

14 J. G. Zabolitzky, Phys. Rev. B 22, 2353 (1980). 

For a recent review see E. Krotscheck and M. Saarela, Phys. Rep. 232, 1 (1993). 

16 T. Chakraborty, Phys. Rev. B 25, 3177 and 26, 6131 (1982); T. Chakraborty. A. Kallio, L. J. Lantto, and P. Pietilainen, 
Phys. Rev. B 27, 3061 (1983). 

17 M. Levy and J. P. Perdew, Phys. Rev. A 32, 2010 (1985). 

18 M. L. Chiofalo, S. Conti, and M. P. Tosi, Mod. Phys. Lett. B 8, 1207 (1994). 

19 M. W. C. Dharma-wardana and F. Perrot, Phys. Rev. Lett. 84, 959 (2000); F. Perrot and M. W. C. Dharma-wardana, Phys. 
Rev. B 62, 16536 (2000). 

20 A. Kallio and J. Piilo, Phys. Rev. Lett. 77, 4237 (1996). 

21 P. Gori-Giorgi, F. Sacchetti, G. B. Bachelet, Phys. Rev. B 61, 7353 (2000). 

22 J. C. Kimball, Phys. Rev. A 7, 1648 (1973). 



8 



D. M. Ceperley, private communication. 

K. S. Singwi, M. P. Tosi, R. H. Land, and A. Sjolander, Phys. Rev. 176, 589 (1968); see also T. K. Ng and K. S. Singwi, 
Phys. Rev. B 35, 6683 (1987). 

S. Tanaka and S. Ichimaru, Phys. Rev. B 39, 1036 (1989); S. Ichimaru, Phys. Rev. Lett. 84, 1842 (2000). 
B. Davoudi, M. Polini, and M. P. Tosi, Solid State Commun. 124, 335 (2002). 



TABLE I: Energy of the 3D EG in Ryd/electron. CA/QMC from Ceperley and Alderyi OHB/QMC from Ortiz et al. 
ZLC/QMC from Zong et al.& STLS from Singwi et and TI from Tanaka and Ichimaru^ 



r s Various calculations £ = 0.0 


C = 0.333 


C = 0.667 


C = 1.0 


1 CA/QMC 


1.174 








OHB/QMC 


1.181 






2.294 


Present work 


1.17810 


1.29786 


1.66346 


2.29753 


STLS 


1.1704 


1.2885 


1.6505 


2.2849 


TI 


1.167 






2.281 


5 CA/QMC 


-0.1512 






-0.1214 


Present work 


-0.14974 


-0.14616 


-0.13557 


-0.11846 


STLS 


-0.1511 


-0.1483 


-0.1393 


-0.1235 


TI 


-0.1544 






-0.1267 


10 CA/QMC 


-0.107 






-0.101 


OHB/QMC 


-0.106 






-0.101 


Present work 


-0.10562 


-0.10488 


-0.10277 


-0.09957 


STLS 


-0.1058 


-0.1055 


-0.1043 


-0.1020 


TI 


-0.1084 






-0.1046 


20 CA/QMC 


-0.063 






-0.0625 


OHB/OMC 


-0.063 






-0.0625 


PrpqpTit, work 


-0 06265 


-0.06250 


-0.06210 


-0.06153 




0fi9^ 
-u.uuzo 






-0.0623 


TI 


-0.0642 






-0.0642 


40 ZLC/QMC 


-0.03523748(60) 


-0.03523295(67) 


-0.03520539(67) 


-0.03513483(72) 


Present work 


-0.03470 


-0.03467 


-0.03459 


-0.03450 


STLS 


-0.0342 






-0.0345 


50 ZLC/QMC 


-0.02889900(62) 


-0.02889962(68) 


-0.02888835(62) 


-0.02884983(81) 


OHB/QMC 


-0.029 






-0.0288 


Present work 


-0.02845 


-0.02844 


-0.02839 


-0.02834 


STLS 


-0.0280 






-0.0282 


TI 


-0.0290 






-0.0293 



0.8 



0.6 



0.4 



0.2 



1 1 
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■ - i i 


Present work 
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FIG. 1: The pair distribution function g(r) in a paramagnetic 3D EG at r s — 1, 5, 10 and 20 (from top to bottom at low r), 
a function of r/(r s as). The results of the present work (full lines) are compared with QMC data of Ortiz et al£ (dots). 
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FIG. 2: The parallel-spin pair distribution function gn(r) in a paramagnetic 3D EG at r s = 1,5,10 and 20, as a function 
of r/(r s os). The results of the present work (full lines) are compared with QMC data of Ortiz et al. 4 (dots). The curves at 
r s — 5, 10 and 20 have been shifted upwards for clarity by 0.4, 0.8 and 1.2, respectively. 
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Present work 
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FIG. 3: The antiparallel-spin pair distribution function <?u(r) in a paramagnetic 3-D EG at r s — 1,5, 10 and 20, as a function 
of r/(r s as). The results of the present work (full lines) are compared with QMC data of Ortiz et al. 4 (dots). The curves at 
r a = 5, 10 and 20 have been shifted upwards by 0.4, 0.8 and 1.2, respectively. 
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FIG. 4: The average g(r) and spin-resolved pair functions in a paramagnetic 3D EG at r a = 10, as a functions of r/(r s Ofl). 
The results of the present work (full lines) are compared with our previous Hartree results- (dashes) and with QMC data of 
Ortiz et al, 4 ' (dots). The results for <?n(r) have been shifted upwards by 0.4 and that for <7fj.(r) by 0.8. 
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Present work 
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FIG. 5: The pair distribution function g(r) in a paramagnetic 3D EG at r s = 40 and 100, as a function of r/(r 3 a,B)- The 
results of the present work (full lines) are compared with QMC data of Ortiz et all— (dots) and of Zong et aZi- (empty boxes) . 
The curves at r s = 100 have been shifted upwards by 0.4. 
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FIG. 6: The pair distribution function g(r) in the fully spin-polarized 3-D EG at r s = 1,4,40 and 100, as a function of r/(r s as). 
The results of the present work (full lines) are compared with QMC data of Ortiz et alm> (dots) and of Zong et oZi- (empty 
boxes). The curves at r s — 4,40 and 100 have been shifted upwards by 0.4, 0.8 and 1.2, respectively. 




FIG. 7: Left Panel (momentum space): the spin-dependent effective interaction nV^f (q) (in units of e 2 /ob) in a 3D EG at 
r s = 10 and 20, as a function of qr s a,B- Right Panel (real space): the spin-dependent effective interaction V^g (r) (in units of 
h 2 kp/m) in a 3D EG at r s = 10 and 20, as a function of r/(r s ae). 
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FIG. 8: The Pauli potential v% a (r) (in units of h 2 k 2 F /m), as a function of r/(r s as). With this choice of units the Pauli 
potential has no explicit dependence on r s . 



